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Approximate  emiscivity  calculations  for  CO,  have  been  carried 
out,  as  a function  of  optical  density,  at  300  and  at  6t>0°K.  The  calculations 
involve  the  assumption  that  the  rotational  lines  overlap  extensively.  This 
condition  appears  to  be  satisfied  at  total  pressures  above  about  1 atmos. 
Comparison  of  the  values  calculated  from  spectroscopic  data  with  the 
emissivities  tabulated  by  fidttel  and  his  collaborators  shows  satisfactory 
agreement.  The  analysis  presented  in  this  manuscript  emphasises  the 
fact  that  it  is  possible  to  obtain  reasonable  estimates  for  the  engineering 
emissivity  without  performing  extensive  analytical  work,  provided  the 
physical  principles  are  understood  and  the  needed  spectroscopic  data  are 
available. 


♦ Supported  by  the  ONR  under  Contract  Nonr- 220(03),  NR  015  210. 
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I.  INTRODUCTION 

In  a series  of  recently  published  papers  we  have  attempted  to 
estimate  the  engineering  emissivities  of  diatomic  molecules  by  utilising 
basic  spectroscopic  constants.  Tractable  methods  of  calculation  have 
been  developed  for  (a)  complete  overlapping  between  rotational  lines 1 and 
for  (b)  separated  rotational  lines.  ^ Although  accurate  theoretical 
calculations  of  gas  emissivities  from  spectroscopic  data  for  polyatomic 
molecules  involve  formidable  computational  difficulties,  it  is  to  be 
expected  that  approximate  calculations  can  be  carried  out  with  relatively 
little  effort  by  utilising  approximations  similar  to  those  introduced  for 
the  study  of  diatomic  molecules . * ' ^ For  CO^  it  has  been  found  that  the 
total  absorptivity  at  room  temperature  is  substantially  independent  of 

3 

pressure  at  pressures  exceeding  about  1 atmos  , thus  indicating 
extensive  overlapping  between  rotational  lines.  Accordingly,  it  is  not 
unreasonable  to  base  the  present  preliminary  calculations  on  the 
assumption  that  extensive  overlapping  between  rotational  lines  does 
occur.  The  results  would  be  expected  to  apply,  for  example,  at  room 
temperature  for  total  pressures  in  excess  of  1 atmos.  and  at  3000°K  at 

3 

pressures  in  excess  of  11  atmos.  Since  it  appears  quite  likely  that  the 


(a)  S.  8.  Penner,  J.  Appl.  Phye.  21,  685  (1950); 

(b)  J.  Appl.  Mech.  |£,  53(1951);  (c)  S.  S.  Penner  and  D.  Weber, 

J.  Appl.  Phys.  22,  1164(1951). 

(a)  8.  8.  Penner,  M.  H.  Ostrander,  and  H.  8.  Tsien,  J.  Appl.  Phys. 
23,  256  (1952);  (b)  8.  S.  Penner,  J.  Appl.  Phys.  ^3,  825  (1952). 

* R.  J.  Holm,  D.  Weber,  and  8.  8.  Penner,  J.  Appl.  Phys.  23,  1283 
(1952).  — 


actual  r*nf«  of  validity  of  the  results  covers  a Larger  range  of  total 
pressures,  we  shall  not  hesitate  to  show  that  calculated  exnissivities 
agree  reasonably  well  with  empirical  data  obtained*  at  1 atmos.  at  a 
temperature#  of  600°K  as  well  as  at  300 °K.  Perhaps  the  most  important 
conclusion  which  can  be  derived  from  the  present  analysis  is  the 
statement  that  the  analytical  labor  involved  in  making  approximate 
emissivity  calculations  for  a polyatomic  molecule  with  overlapping 
rotational  linos  is  trivial  provided  the  important  physical  principles  are 
understood  and  the  needed  spectroscopic  data  are.  available.  It  is  only 
fair  to  note,  however,  that  all  of  title  needed  integrated  intensity 
measurements  have  not  yet  been  performed  for  CO^  in  spite  of  the 
fact  that  several  papers  dealing  with  intensity  measurements  have  been 
published  recently.5,6, 7 

n.  basic  relations 

Theoretical  considerations  of  intensities  for  various  transitions 
of  the  CO2  molecule  were  summarised  by  Dennison  a number  of  years 
ago. 6 We  shall  reproduce  here  the  parts  of  the  analysis  which  are 
useful  for  making  approximate  emissivity  estimates  on  CO^. 


4 W.  H.  McAdams,  Heat  Transmission.  McGraw-Hill  Book  Co. . 

New  York  1942,  Chapter  111  by  H.  C . Hottel. 

5 A.  M.  Thorndike,  J.  Chem.  Phys.  15,  868  (1947). 

6 D.  F.  Eggers,  Jr.,  and  B.  L.  Crawford,  Jr.,  J.  Chem.  Phve.  19, 
1554(1951). 

7 D.  Weber,  R,  J.  Holm,  and  S.  S.  Penner,  J.  Chem.  Phys.  20, 

1820(1952).  “ 

8 D.  M.  Dennison,  Rev.  Mod.  Phys.  3,  280  (1931). 
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FrMjaMciti  (O)  and  wav*  numbers  (O)  corresponding  to 
transitions  between  fixed  energy  levels  are  given  by  the  Bohr  frequency 
condition 

O • c U * (Atr'y  + dw'R)/h  <1) 

t « 

where  AW  y and  AW  ^ represent,  respectively,  the  changes  in 
vibrational  and  rotational  energy  corresponding  to  the  frequency  U, 
h is  Planck's  constant,  and  c represents  the  velocity  of  light.  According 
to  the  results  of  Dennison  and  Adel,  the  rotational  (W  and 

9 

vib rational  (W  y)  energies  for  C02  are  given,  respectively,  by  the 
expressions 

w'R  * hc(J2  + j -lZ)  [0.  3925  - 0.  00058  (nj  + -J) 

+ 0.  00045  <n2  ♦ 1)  • 0.  00307  <n3  ♦ ■£>]  +1.7  - 1)  (2) 

and 

w’y  <•  he  [l351.  2 (nt  + -J)  + 672.  2 (n2  + 1)  ♦ 2396.  4 (n3  + 4) 

- 0.  3 (nj  + ^j)2  - 1.  3 (n2  + l)2  - 12.  5 (n3  + -J)2 

0 

+ 5.7  (nj  +^)  (»2  ♦ 1)  - 21.9  (nj  +4'  (»3  ♦ \) 

- 11.0  (n2+  i)(n3+,2)]  . (3) 


D.  M.  Dennison,  Rev.  Mod.  Phys.  12,  175  (1940). 

Slight  corrections  to  the  listed  numerical  values  have  been  noted  by 
W.  S.  Benedict  and  his  collaborators.  For  the  present  purposes 
these  corrections  are  negligibly  small.  Recent  work  on  frequencies 
of  lines  belonging  to  CO,  is  described  in  the  following  papers: 

W.  S.  Benedict,  S.  C.  Herman,  and  S.  Silverman,  J.  Chem.  Phys. 
19,  1325  (1951);  J.  H.  Taylor,  W.  S.  Benedict,  and  J.  Strong, 
Progress  Report  on  "Infrared  Spectra  of  H,0  and  CO,  at  500°C", 
Contract  Nonr.  248-01,  The  Johns  Hopkins  University,  March  1952. 


Here  , a^,  end  ere  the  vibrational  quantum  numbers  associated 
with  the  fundamental  vibration  frequencies  \?j,  V>2,  and  Vj.  respect- 
ively. The  quantum  number  X measures  the  angular  rotation,  in  units 
of  h/2v,  which  is  associated  with  the  degenerate  (bending)  \?2~  vibration. 
The  rotational  energy  levels  are  identified  by  the  rotational  quantum 
number  J. 


Q / 

The  selection  rules  for  the  transitions  n^  n? 


. I' 


n , n2  n3 


may  be  summarised  as  follows.  For  the  perpendicular  bands,  An^  odd, 
An^  even,  A X ■ + 1.  For  the  parallel  bands,  An^  even,  An^  odd, 

A^  «0.  The  rotational  selection  rules  are  Aj  = 4 1,  0. 

8 11 
The  amplitude  factors  for  given  rotational  transitions  are 


(Aj  if  *tt mr 

(4) 

(A]-i  t)  ' Wn)  * (Aj  s) 

(4a) 

(4b) 

W 

(4c) 

<4d) 

i i ) 8 rferi)  • 

<4e) 

vr*  Aj  •**$&&*•**  * 

(4f) 

The  amplitude  factors  given  in  Eqs.  (4)  to  (4g)  are  4 times  as  large 
as  those  listed  by  Dennison**  because  we  are  using  a Fourier  series 
in  time  of  the  form  cos(2ri/t)  rather  than  exp  (2«i<?t). 


v 
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(*i • ' » • Mb,  ■ 


Hi) 


b Eqi.  (4)  to  (4g)  tho  convention  h&i  boon  adopted  that  tho  amplitude 

M * M ' 

factor  corresponding  to  the  transition  J,X--»  j , / has  been  written  as 

M i)- 

The  Integrated  abeorption  for  a (Ivan  transition  from  the  lower 

i 1 • jC  « • 

energy  level  njU^  a^l  j to  the  upper  energy  level  n j n ^ n y,  J 

JL  • t i • g 

is  designated  as  8 m Sfnj^  n^;  j — *n  j n 2 n 3;  J ) and  is  given  by 


8inxn 

the  approximate  relation 


<8v3^NT/5hcQ'v  q'r)|  exp  - [w'y^.  n^.  n3)  ♦ w'R  (),  J )]  /Mr} 

hv>, 

*'  • • 


*«i 


fl2  [l  - exp  (-hv>/kX)l 

\ i i ) 

L J 

iS) 


n 1 n2  n y J ) is  the  Bohr  frequency 


weight  of  the  upper  state  with  Sj«y i * IJ  ♦ 1 for  t ■ 0 and  gj , ,,,  « 2(2j  + 1) 

#'  / i » 

for  Jb  ^ Oi  f gp  (n j n^*  n,  —a  n , n7na  4)  is  a factor  which  must 


In  *q.  <5)  »s  \>(n|  n^  n^  J 
corresponding  to  the  indicated  change  in  the  quantum  numbers;*  NT  « total 
number  of  molecules  per  unit  volume  per  unit  pressure;  %yj^\  * statistical 

*JV 

t 

“3  n 1 “2  “3 

be  determined  empirically  and  corresponds  to  the  matrix  element  of  the 
electric  moment  in  the  molecule  associated  with  the  indicated  change  in 

(vibrational)  quantum  numbers;*  Q R TZ  gj^  exp  W R(J,</)/kTj  » 

complete  rotational  partition  function; 

o'vE  n z=  exp  - [W  y(np  n^,  n3)/kX  j * complete  vibrational 
nl  nZ  n3 
partition  function. 

* Although  we  shall  not  indicate  explicitely  the  changes  in  the  (vibrational) 
quantum  numbers  involved,  both  v’ and  0 do,  of  course,  vary  from  one 
vibration- rotation  band  to  another. 


( > 

* 


r 


* 


The  amplitude  factors  given  in  Eqs.  (4)  to  (4g)  are  based  on  the 
assumption  that  vibration -rotation  interactions  have  a negligibly  small 
influence  on  amplitude  factors.  For  this  reason,  the  use  of  Eq.  (5) 
involves  somewhat  cruder  calculations  than  were  employed  in  the 
calculations  on  diatomic  molecules  with  non-overlapping  rotational  linos.  * 
For  the  purposes  of  approximate  radiant  heat  transfer  calculations 
it  is  convenient  to  use  the  integrated  absorption  coefficient  for  a given 
band.  Lot 


J • • v 
«i“l“2  V nln2  n 


-£  E 

1 i y 


• » J 
1 J2 


n ,n  , n J ).  (6) 


For  the  parallel  bands  di  * 0 and Rj'/'  ^ 
according  to  the  Burger  and  Dorgelo  summation  rules,  which  may  be 
verified  by  use  of  Eqs.  (4)  to  (4g)  by  utilising  the  selection  rules.  Hence 
obtain  from  Eqs.  (5)  and  (6)  the  useful  result 


• i • • ( - 

s^nj^'^Uj  —v  n jn  2 n j)  * (8wSp2  *NT/3hcQy)gJ  exp  -IWyinj.n^nj.JPj/kT]  > 

* 


t[l  - exp(-hWkT)|  with  gg  * 1 for * 0 and  2 
for  0.  (7) 


where  we  have  set  Wy*  (n^n^n^)  + WR*  (J.i)  * Wy<n|.n2.n),/)  + WR(J)» 

Qv  * ,xp  <-wv/«).qv-qr-  . qvqr  With  Qr  -E  «j  «p[-WU>/m]. 

and  gj^£  ■ gj  x g j . If  the  assumption  is  not  made  that 


• * 


>)•  \>'ar  V (®|n2  ^ 0 -♦  n jtt  ^ a j|  0)  e alight  correction  1c 

obtained  to  Eq.  (T). l2* 13 

Tor  the  perpendicular  bands  A-/ + 0 and  t.  - - gj,^,  <xj'  ^}2  « g^ 

j'*/  1 

IT  <Aj  * -j  «J/’  5JIZ  S^,  (a}  V -fj^.  Hence 


^“iV^  -*  n\  * ln's)  • 


(4*V  i>'NT/3bcQv)  g^ 


|e*p  - [Wy<nl,n2.B3,|)/kTj||i  -exp<-h\>/kT)J 


where  g^  • 1 for  0 and  g^  * 2 lorjtit  0,  (8) 


*<miV**3 


V3)  « (4»V  ^Nj/ShcQyJ  g^ 


x^axp  -Jw^nj.ng.nj./j/kTjj  Jl  - exp(-hV/kTf] 
where  g^  * 1 for  A*  0 and  g^  * 2 for  jt+  0.  (9) 


In  using  Kqs.  (8)  and  (9)  it  should  be  noted  particularly  that  the  quantities 
f are  generally  different  for  the  transitions  involving  ^ ~+u£ ♦ 1 and 


B.  L.  Crawford,  Jr.,  and  K.  L.  Dina  more,  J.  Chem.  Phys.  18.  983, 
1682  (1950). 

Bggers  and  Crawford^  have  shown  how  to  relate  for  CO,  the  matrix 
coefficients  to  the  coefficients  of  Taylor  series  expansions  of  the 
potential  energy  and  dipole  moment  about  their  respective  equilibrium 
values. 
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I.  Tht  ratios  of  th«  p*  con  bo  calculated  theoretically  to  tho 
harmonic  oscillator  approximation  for  harmonic  bands  (so#  Appendix  I 
for  details). 

In  the  same  manner  as  for  emissivity  calculations  on  diatomic 
molecules  with  overlapping  rotational  lines,  * we  shall  find  Eq.  (5) 
useful  for  the  determination  of  “effective  band- widths"  whereas  Eqs.  (7) 
to  (9)  lead  directly  to  "average  absorption  coefficients"  for  vibration- 
rotation  bands.  For  the  calculation  of  S it  may  be  convenient  to 
combine  Eqs.  (4)  to  (4g),  (5),  and  (7)  to  (9).  For  example,  for  the 
positive  branch  of  the  - fundamental  of  CO^.  we  obtain 

S **  JoXOlV  j -*00°0;  J-l)/  0>  (O^OjO  ► 00°0;0)]  x <QRf 1 

x —a  00*0)  x (j  + 1)  [exp  - <WR/kT)j  (5a) 

where  GJ  denotes  a wave  number. 

m.  approximate  emissiyitt  calculations  at  3oo°k 

The  positions  and  approximate  intensities  of  the  stronger  vibration- 
rotation  bands  ol  CO^  at  3C0°K  are  well-known.  ^ A first  approx- 
imation to  the  effective  band-width  is  obtained  by  utilising  relations  such 
as  Eq.  (5a)  and  defining  the  effective  band-width  as  the  wave  number 
range  for  which  S exceeds  10'  of  its  maximum  value.  Utilising  this 
definition  of  band- width,  which  we  shall  refine  presently,  a summary  of 
vibration- rotation  bands  for  CO2  at  300°K  has  been  drawn  in  Figs,  i to  3. 

• 

**  G.  Hersberg,  Infrared  end  Raman  Spectra  of  Polyatomic  Molecules, 

D.  Yan  Nostrand  Co. , View  York  1951,  Table  56  on  p.  Z 74.  HersEerg's 
designations  of  band  intensities  as  vs  (very  strong),  s (strong) 
m (medium),  w (weak)^and  vw  (very  weak)  has  been  used  in  Figs.  1 to  3 
to  identify  the  approximate  strengths  of  bands  at  room  temperature. 


Table  I.  Integrated  intensities  for  vibration  - rotation  bands  of  C02. 7 


Band  Center 

(cm-1) 

Transition 

2 1 
0C(cm”  -atmoi  ) 

5109 

4984 

4861 

3716 

3609 

2349 


1933 

1886 


961 

1064 


01 10 
0I10 


— * 20°0  ^ 
— 12*0 
— * 11*0  J 


01l0 

oilo 

01*0 

00°0 

02°0 

01l0 

10°0 

02°0 


10°0 

0220 

01*0 

03*0 

02°0 

00°  1 

00°1 


0.426 

1.01 

0.272 

42.3 

28.5 

2706 

. 147 
.083 
171.5 


0.0219* 

0.0532* 


* Unpublished  dsta  obtained  by  O.  Web  jr 
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7 

as  they  are  available.  For  the  strongest  bands  (O2  ~ *®d\^  - funda- 

13  16 

mentals)  the  contributions  from  the  isotopic  molecule  C are 

also  shown.  This  chemical  species  constitutes  about  1. 1 °/o  of  all  of 
the  CO^  present  and  makes  contributions  to  the  total  emisaivity  which 
are  not  negligibly  small  for  moderate  and  large  optical  densities  of 
emitter.  Also  shown  in  Figs.  1 to  3 are  the  spectral  intensities  of 
radiation  emitted  by  a blackbody  at  300,  600  and  1500°K,  again  utilizing 
arbitrary  scales  along  the  ordinates,  which  are,  however,  consistent  for 
any  one  of  the  given  temperatures. 

Reference  to  Fig.  1 shows  that  the  principal  contribution  to  the 
total  radiant  heat  transfer  at  room  temperature  must  occur  from  the 
spectral  region  around  the  ^ ~ fundamental.  Thus  it  is  readily  shown* 
that  the  "limiting  emissivity”  for  CO2  at  300°K  is  in  the  neighborhood 
of  0.  4 with  more  than  90  °/o  of  the  total  contributions  made  by  the  wave 
number  region  extending  from  551  to  830  cm  *. 

It  is  physically  obvious  that  the  engineering  emissivity  E at  300°K 
is  determined  by  the  01*0  — ► 00°0  transition  in  the  limit  of  zero  optical 
density  of  the  emitter.  As  was  pointed  out  in  an  earlier  publication,  * 
the  emissivity  as  a function  of  optical  density  (partial  pressure  p of 
emitter  z optical  pathlength  L)  for  very  small  values  of  p L can  be 
computed  quite  accurately  either  by  utilizing  a treatment  applicable  for 
non -overlapping  lines  ^ or  by  using  average  absorption  coefficients.  The 
emissivity  of  the  band  with  center  at  667  cm'*,  as  calculated 

from  Eq.  (30)  of  Ref.  2a,  is  plotted  as  a function  of  pL  in  Fig.  4 and 
labelled  "approximation  for  non-overlapping  lines".  The  corresponding 


ll  < 


•i2- 
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quantity  obUlstd  by  using  the  average  absorption  coefficient 
PX  * (17i  cm"2  atxnos”  *)/( 106  cm”  *)  over  the  effective  band-width 
extending  from  600  to  706  cm" * give*  results  which  are  practically 
identical  with  those  obtained  for  non-overlapping  rotational  lines  as  pL 
is  reduced  below  about  0. 02  ft-atmos . As  pL  is  increased  Eq.  (30)  of 
Ref.  (2a)  rapidly  fails  to  apply  and  yields  excessively  large  values  not 
only  for  bat  also  for  the  total  engineering  emissivity  E.  Need- 

less to  say,  for  pL  <0.02  ft-atmos,  the  numerical  value  of  E^7  is 
to  be  identified  with  the  value  of  E. 

As  the  optical  density  is  increased  to  moderate  values  of  p L, 

it  bee  omen  necessary  first  to  include  the  01*0  » 00°  0 transition  for 

C**  which  is  assumed  to  represent  1.  1 °/o  of  the  total  COj. 

and  for  which  the  integrated  intensity  is  roughly  the  same  as  for  the 
12  16 

C <>2  species.  In  treatments  using  average  absorption  coefficients 
it  is,  of  course,  necessary  to  add  absorption  coefficients  rather  than 
emissivities  in  regions  in  which  partial  overlapping  between  effective 
band-widths  occurs.  Finally^ ths  contributions  made  by  the  10° 0 — ► 01*0 
and  02°0  ■ -■«  01  *0  bands  must  be  included;  both  of  these  are  designated 
by  Hersberg**  as  being  of  "medium"  intensity  whereas  the  03*0  — ► 02°0 

and  11*0  — — ♦ 0220  bands  are  labelled  "weak"  and  the  11*0  *02°0 

band  is  said  to  be  "very  weak" . 

Strictly  speakiug) room  temperature  emissivity  calculations  for 
CC>2  cannot  be  carried  farther  because  the  integrated  intensities  for  the 
weaker  bands  are  not  available.  However,  it  is  easily  shown  that 
results  in  fairly  good  agreement  with  Hottel’s  data.*  are  obtained  if  reason- 


-n- 


abl#  estimates  an  mad*  for  th*  *trong**t  of  the  band*  which  has  not 
yet  been  included  in  determining  S. 

7 

It  ie  known  from  rough  unpublished  measurements  that  the 
integrated  intensity  ^721  for  the  i0°0  — — *0iJ0  band  is  of  the 
order  of  1 cm" * atmos”*.  Emissivity  calculations  treating  ^721  *• 
a variable  parameter  between  0 and  10  cm”*-atmoe”*  are  shown  in 

Fig.  4.  The  corves  labelled  ^721  * °*  ^721  * 1111(1  °^72l  0 10 

with  0(  ^ *171  correspond  to  the  calculated  values  of  E for  the 
various  assumed  values  of  pL.  In  making  these  calculations  tha  concept 
of  th*  effective  band- width  was  refined  in  so  far  as  the  effective  band-width 
for  th*  fundamental  was  set  equal  to  the  wave  number  range  for 

which  SpL/0.7  >0. 1.  This  definition  of  band-width  for  strong  bands 
has  the  important  advantage  of  making  the  band- width  a weak  function 
of  pi*,  in  agreement  with  empirical  measurements.  The  contributions 
to  E made  by  the  00°  1 —*■  10°0  and  00°1  — ► 02°  0 bands  are  very 
small  but  have  also  been  included.  To  these  bands  Weber  assigns  the 
following  values  at  300°K:  0219  cm”*  atmos”*, 

0(  • 0.0532  cm”*atmos”1.  

In  addition  to  the  calculated  emissivities,  the  results  of  empirical 
measurements*  are  also  shown  in  Fig.  4.  For  pL  ^ 0. 06  ft-atmos. 
the  nempiricalH  data  are  really  extrapolated  results.  Reference  to  Fig.  4 
shows  fair  agreement. between  calculated  and  obeerved  values  of  E for 
the  reasonable  values  *667  « 171  and  0(  J21  * 1 czn  ^ataos  \ At 
very  small  values  of  p 1*  the  calculated  emissivities*  are  larger  than 

* The  dotted  curve  in  Fig.  2 corresponds  to  / / - = 29  and 

_2i  * 1 cm~2atmos~l.  Although  this  curve  represents  a good  fit 
of  tht empirical  data,  the  result  cannot  be  considered  to  be  - 
significant  since  *667  is  known  to  be  much  larger  than  29  cm”  atmos 


1 


t 


the  oUtpolatad  MB^vlcal  data  bat  mast  be  coaiid«r«d  to  be  mott 
rritobU  iliac*  thiy  dipnd  only  on  the  numerical  mhu  of  <K  , which 

U known  with  fair  accuracy.  * Oar  ability  to  calculate  E at  vary 
largo  vatnao  of  p L ii  lixnitad  by  tha  lack  of  adequate  intensity  data  for 
tha  tramaitiooa  02°0— *0Il0,  11*0  — *0*20.  and  11X0  — *02°0.  How. 
ever,  it  is  clear  that  ai  pL  ii  incraaaad  sufficiently  abova  9 ft-atmos, 

E moat  approach  about  0. 4 and  will  than  increase  only  vary  slowly  as 
aacaadlagly  largo  vaboes  of  pL  ara  ebOinad. 

Tha  discussion  of  amissivity  calculations  on  C02  at  300°K  as  a 
function  of  pL  is  Mostly  of  academic  interest.  However,  it  is 
instructive  in  so  far  as  it  exampliflee  tha  inherent  simplicity  of  tha 
analysis  if  adequate  spectroscopic  data  are  available. 

C IV.  APPEOQOhiATS  EMISSIV1TT  CALCULATIONS  AT  600°K 

It  is  evident  from  Fig.  2 that  tha  major  contributions  to  radiant 
heat  transfer  at  600°K  ara  made  by  vibration- rotation  bands  in  tha 
spectral  region  between  >90  and  2400  cm*1.  As  pL  goes  to  sero  it  is 
expected  that  the  total  amissivity  E will  become  substantially  equal  to 
the  amissivity  of  the  intense  \?  ^ * fundamental  of  with 

appropriate  corrections  for  contributions  from  the  isotopic  species 
C^O^.  Although  we  are  unable  to  complete  the  amissivity  calculations 
at  600°K  becauss  the  needed  spectroscopic  dsta  are  not  available,  the 
engineering  amissivity  is  estimated  correctly  for  small  values  of  p L. 

This  result  emphasises  the  fact  that  relatively  accurate  amissivity 
calculations  for  polyatomic  molecules  can  be  performed  without  difficulty 
as  soon  as  the  necessary  infrared  intensity  measurements  are  available. 
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For  the  «Mk«r  band*,  i.  •. ; for  all  bands  but  the  'Oj  - fundamental, 

we  shall  define  the  effective  band-width  as  the  wave  number  region  for 

-3 

which  S exceeds  about  10  of  its  maximum  value.  For  the 
VJ  j - fundamental  we  set  the  effective  band-width  equal  to  fixe  wave  number 
range  for  which  (SpL/0.7)  > 0.  i.  * The  long  wave-number  limit  of  the 
Vj  - fundamental  is  set  equal  to  2410  cm’1,  i.e. , 11  cm’1  beyond  the 
band  head.  This  choice  allows  for  fixe  many  tails  of  rotational  lines  which 
"spill"  across  the  band-head  limit.  We  proceed  by  calculating  separately 
the  partial  emissivitiee  for  selected  wave  number  regions. 

(A).  Contributions  of  the  0C°1  -» 10°0  and  00°  1 -»  02°0  Bands 

The  criterion  (S/S^^)  > 10  for  the  weaker  bands  leads  to  band- 
widths  of  fixe  order  of  140  cm'  * at  600°K.  The  integrated  intensity  for  a 
given  vibration-rotation  band  at  the  arbitrary  temperature  T,  divided  by 
the  measured  value  at  300°K,  is  obtained  by  the  use  of  Eqs.  (7),  (8)  or 
(9).  From  these  relations  we  find  for  fixe  integrated  intensity  of  the 
band  whose  ground  state  is  identified  by  the  set  of  quantum  numbers 
' Uj  fixe  result 


♦ Although  the  definitions  of  the  effective  band-widths  are  somewhat 
arbitrary,  it  is  aasily  shown  that  calculated  values  of  E are  quite 
insensitive  to  fixe  chosen  bead- widths.  This  result  is  caused  by 
automatic  choice  of  too  small  absorption  coefficients  for  band -widths 
which  are  too  large,  and  conversely  (Compare  Ref.  la).  Use  of  the 
quantity  SpL/0.7  >0.1  assures  inclusion  within  the  band -width  of 
the  K,  - fundamental  of  all  of  the  rotational  lines  with  average 
emissmties  in  excess  of  0. 1.  Substantially  the  same  numerical 
values  for  the  effective  band-widths  of  the  ^-fundamental  are 
obtained  if  fixe  lines  with  SpL/0.7  > 10“ *n  are  included,  where  n is 

a number  which  does  not  differ  from  unity  by  more  than  a factor  of 
three  or  four. 
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0MT“K)/*(<300oK)  • (300/T)  £ Qy  (300°K)  exp  [wy  (00°0)/300  k]^ 

* ^Qy  <T°K)  * cap  [wy  (00°0)/k7  jj  1 
x |exp  - [wY  (njn/aj)  - ^00°0)][(kTf 1 - (300  k)'l]j 
x £l  - exp  (-h\^kT)j  |l  - ®*P  (-k  9/300k)J  (10) 


where 


• 1 
-2 


Qy(T)  exp  Wy(00°0)/kTj  ^ Jl  - exp  (.1388  kc/kT)| 

x [l  - exp  (-  667  hc/kT)] 
x[l-exp  (-2349hc/kT)]  . 


Cll) 


The  om  of  Eqi.  (10)  aad  (11)  load*  to  the  following  result*  for  the 
00°1  — *»  10°0  and  00°l — *02°0  bands,  respectively:  <*961<*00°K)  — 

9.29  0<961  (300°K)  « 0.204  cm’2  atmos'1;  U 1(J64(600°K)«  7.22  * 1064<300°K) 
■ 0.  385  cm~2atmos~*.  Hence  the  total  contribution  to  the  emissivity 
obtained  from  the  two  bands  under  discussion  is 


E96l  * E1064  “ 0,052 

- exp  (1.40  x 10~3  pL)  J 

4 0.017  | 

1 - exp  (4. 15  x 10' 3 pL)J 

+ 0. 054 

[l  - exp  (2.75  x 10’S  pL)] 

where  pL  is  expressed  in  cm -atmos. 


(B).  CoBtribattw  Made  by  th»  Banda  Bftwtw  1800  and  2210  cm' 1 

It  U known  from  the  recent  work  of  Taylor,  Benedict,  and  Strong15 
that  an  appreciable  number  of  vibration- rotation  bands,  in  addition  to  the 
bands  shown  in  rig.  1 in  this  region,  are  observable  at  500°C.  * Since 
all  of  these  bands  involve  transitions  to  excited  energy  levels,  their 
integrated  intensity?  would  be  expected  to  increase  as  the  temperature  ia 
raised.  Rough  absolute  intensity  estimates  for  several  bands  in  the 
5 micron  region  are  given  in  Table  Q at  300  and  600°K.  The  relative 
intensity  estimates  were  obtained  by  using  the  greatly  oversimplified 
harmonic  oscillator  approximation  described  in  Appendix  1.  Thus  the 
contributions  of  both  mechanical  and  of  electrical  enharmonic! ties  to  die 
intensity  have  been  neglected. 

As  is  evident  from  the  experimental  data  shown  ia  Reference  IS, 
a large  number  of  vibration- rotation  bands  remain  for  which  even  rough 
intensity  data  are  not  available.  For  this  reason  no  result  can  be  given 
for  the  contributions  mads  to  the  engineering  emieeivity  by  die  vibration- 
rotation  bands  between  1800  and  2100  cm~*. 


J.  H.  Taylor,  W.  S.  Benedict,  and  J.  Strong,  Pro greee  Report  on 
"Infrared  Spectra  of  H,0  and  CO,  at  500°C",  Contract  Nonr  248-01, 
The  Johns  Hopkins  University,  Baltimore,  hid. , March  1952. 

* Data  of  ths  type  presented  in  Ref.  15  are  valuable  aide  in  performing 
emieeivity  calculations.  For  example,  the  spectroscopic  Manumission 
records  indicate,  at  least  qualitatively,  the  vibration- rotation  bands 
which  mgst  be  included,  at  temperatures  up  to  773°  K. 


TabU  11.  Absolute  intensities  at  300  and  at  600°K  for  several  vibration- 
rotation  bands  in  the  S micron  region,  based  on  intensity 
calculations  using  harmonic  oscillator  approximations  (Cf. 
Appendix  1.  ).  The  value  0. 147  cm~*-atmos  for  the  sum  of 
the  bands  centered  at  2137,  2094,  and  2077  cm~l  served  as 
standard  to  fix  the  absolute  intensity  scale. 


Transition 

0(at  300°K  (cm"2-atmos”  *) 

o(at  600°K  (cm”2-atmo«" 

oo°o  — ► li'o 

0.12 

0.040 

01X0— ► 1220 

0.020 

0.  032 

01l0— * 20°0 

0.0045 

0. 0075 

01  *0  — *-12°0 

0.00048 

0.  00079 

0220  — * 1330 

0.0012 

0.0099 

0330  — * 1440 

0.000038 

0.0020 

0220— - •2110 

0.00019 

0.0015 

IO°o — •>  21*0 

0.00029 

0.0026 

lll0  — ►so0© 

0.0000053 

0.00026 

(C).  Contributions  from  the  V 3 * Fundamental 

From  the  measured  integrated  intensity^  for  the  V ^-fundamental 
at  300°K  and  by  use  of  Eqs.  (10)  and  (11)  we  find  °^2349  * *060  cm  2itmos*’ 
at  600°K.  Of  the  total  integrated  intensity  per  (cm"  S*(cm*  *-atmos”*  of  CC^) 
we  associate  98.9  °/o  with  Cl20*6  and  1. 1 °/o  with  C13  O^6.  For  the 
intense  V ^-fundamental  the  effective  band- width  corresponds  to  the 
wave  number  range  £4)2349  *or  ^bich  (SpL/0.  7)  >0. 1 with  the  band- 
width of  the  isotopic  species  determined  similarly,  hi  general,  the 


Cl20*6  - band  U black  before  the  isotopic  band  makes  appreciable 
contributions.  The  emissivities  at  600°K,  &2349*  c*lculated  by  the 
use  of  average  absorption  coefficients  P « tf/ACJ,  are  summarized  in 
Table  IQ  as  a function  of  pL.  Comparison  of  ^£349  ^total)  with  the 


Table  m.  Emissivities  E2349  for  C12<>26  and  Cl3Ol26  at  600°K  as  a 
function  of  optical  density.  * 


pL,  cm -atmos. 

pL,  ft -atmos. 

E2349(C  °2  ' 

*tJ49«cUo2t>  E234,<to“1» 

0. 1 

0.0033 

0.018 

0. 0005 

0.019 

0.5 

0.0164 

0.033 

0.0016 

0.035 

1.0 

0.033 

0.037 

0.0025 

0.040 

5.0 

0. 164 

0.050 

0. 0070 

0.057 

15 

0.492 

0.051 

0.011 

0.062 

50 

1.64 

0.052 

0.011 

0.063 

100 

3.28 

0.053 

0.012 

0.065 

200 

6.56 

0.055 

0.012 

0.067 

* The  tabulated  values  of  E£349(cl3°2^  corr*»P<md  to  values  which 

must  be  added  to  E23^(C^202^ ) in  order  to  obtain  E£349  They 

do  not  include  emissivities  in  the  (black)  region  of  the  9^ -fundamental 

ofC*V‘ 

which  is  overlapped  by  the  Vj. 

fundamental  of 

C*V‘. 

empirical  data*  plotted  in  Fig.  5 shows  that  for  pL<  0. 033  ft.  -atmos, 
^2349  i*  practically  identical  with  the  engineering  emissivity  E. 
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The  result  la  in  accord  with  expectations  and  shows  clearly  that  useful 
emlssirity  data  can  be  obtained  with  practically  no  analytical  work 
once  the  physical  principles  are  understood  clearly. 


(D).  Contribution  of  the  02°  1 — » 00°0  and  10°1  — *-00°0  Bands 


From  the  measured  integrated  intensities  we  find  o(  = 
9.70  em’^stmos’^  for  the  02°  1 — s00°0  -band  and  5 

14,4  cm”2atmos"  * for  the  10°1— *00° 0 - band  at  600°K.  The  effective 
band-widths  extend  from  3539  cm"1  to  3679  cm*1  and  from  3646  cm-1 
3716  cm"1,  respectively.  Using  average  absorption  coefficients  the 
following  results  is  obtained: 


E3609  + 


‘3716 


- 0.  0047  [l  - exp  (-6. 92  x 10" 2 pL)J 

♦ 0.0017  £ i - exp  (-1.72  x 10_1  pL)] 
+ 0.  0017  [l  - exp  (-1.03  x 10" 1 pL)1  . 


(E).  The  Total  ( Engine e) 


Emissivity 


The  total  (engineering)  emissivity  is  obtained  by  adding  the  partial 
emissivities  from  selected  spectral  regions.  Since  estimates  of  partial 
emlssivities  were  not  included  for  all  of  the  important  vibration-rotation 
bands,  no  general  expression  for  the  engineering  emissivity  can  be  given 
although  useful  results  have  been  obtained  for  small  values  of  the  optical 
density  (compare  and  Fig.  5). 


// 
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APPENDIX  I.  RELATIVE  INTENSITY  CALCULATIONS  TOR 
COj  USING  HARMONIC  OSCILLATOR 
APPROXIMATIONS* 


St  is  wtll  known  that  harmonic  bands  with  aon-sero  matrix  com* 
poneste  art  predicted  am  to  tha  harmonic  oscillator  approximation. 
Caicnlationa  of  this  sort  ara  useful  for  rough  r slat! vs  intensity  estimates 
although  mechanical  as  well  as  electrical  anharmonlcity  corrections  are 
neglected.  Representative  relative  intensity  estimates  have  been  given, 
for  exapnple,  by  Benedict**  and  by  Kaplan.  **  For  the  present 
applications  it  will  be  conveniens  to  present  a general  equation  tor 
relative  intensities  of  harmonic  bands  involving  changes  in  the  vibrational 


qpaatam  number  \)j. 


The  desired  result  is  obtained  readily  by  using  the  normalised 

wave  functions  tor  the  isotropic  plane  oscillator  and  an  integral  involving 

e*  17 

associated  La guerre  polynomials  first  obtained  by  Schrodinger.  In 
this  manner  it  is  easily  shown  that 


£ (n1n2'*a3-*.n1n  * Uj)J  ^ [pinjO^-enjl*!^)]2  « Rln^-yn  2 ^ ) • 


4 (^XrfX'rt' 


(A- 1) 


* The  author  is  indebted  to  Drs.  W.  8.  Benedict  and  L.  D.  Kaplan  for 
helpful  correspondence  concerning  intensity  estimates  for  harmonic 
binds  of  CO^. 

16  L.  D.  Kaplan,  J.  Cham.  Phys.  _18,  186  (1950). 

17  E.  Schrodinger,  Ann.  d.  Physik.  80.  483  (1926). 


where  k » (1/2)  (n?-^),  k'  ■ (1/2)  (n'2  -j£)t  p « (1/2)  (J  + 2 ♦ 1).  And 
denotes  a binomial  coefficient. 

By  combining  Kqa.  (7)  to  (9)  with  Eq.  (A-l)  the  following  useful 
result  is  obtained: 


<X <nin2 f °s  nin  2^  ®S^  <®|  °°n3  — ► 1 «3)  s ) 

/ • j,'  r 2 ' jI'  , o 1 

* R(»2  “^n2  n3~*nln2  1^0  nj-^njlB 


:|«sp  - ^Wy  (nj.n^nj,^)  - Wy  (nj,  0.  n^,  0)  j /kT^ 


j\  - exp 


h?(n,  <en 


^ )/kTjj*|l  - exp  £-hV(  0°“*  H/kT^  . (a- 


j . « jl 

Numerical  values  of  Rjn^  -fn  2 ) calculated  from  Eq.  (A-l)  are 

15 

identical  with  Benedict's  estimates,  which  utilised  data  given  by 
Shaffer. 1 


Following  Benedict  we  shall  assume  that  resonance  between  2 V2 
and  Vj  1*  exact,  an  approximation  which  will  not  change  the  order  of 
magnitude  of  the  calculated  results  although  it  is  not  in  accord  with 
observed  infrared  and  Raman  intensities.  By  utilising  Eqs.  (A-l)  and 
(A -2)  relative  intensities  have  been  computed  for  a number  of  bands  in 
the  5 )*  region  at  300  and  at  600°K.  The  results  are  summarised  in 
Table  (A-l). 
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W.  H.  Shaffer,  Rev.  Mod.  Phys.  16,  245  (1944). 


LIST  or  FIGURES 


Fig.  1 . PmMom  of  centers  of  vibration* rotation  bends  observed 

et  room  temperature  between  300  end  2400  cm’1.  The 
designations  vs  (very  strong),  s (strong),  m (medium), 
w (week),  end  vw  (very  week)  ere  those  of  Hersberg.1* 
Also  shown  is  the  intensity  ratio  R^/R^  for  e 
bleckbody  et  300°K. 


Fig.  2.  Positions  of  centers  of  vibration- rote tion  bends  observed 

et  room  temperature  between  600  end  4800  cm’1.  Also 
shown  is  the  intensity  ratio  R ^ /R^  for  e bleckbody 

et  600°K, 


Fig.  3.  Positions  of  centers  of  vibration-rotetion  bends  observed 

et  room  temperature  between  1500  end  12,  000  cm’1.  Aleo 
shown  is  the  intensity  ratio  / R)^  n)aT  for  e bleckbody 

et  15ou°K. 


Fig.  4.  Engineering  emissivity  E as  e function  of  pL  for  C02 

et  300°K.  Integrated  intensities  X ere  expressed  in 
cm’^-ataos’1.  The  results  based  on  empirical 
measurements  (Cf.  reference  4)  are  shown  as  circles. 

Fig,  5.  Engineering  emissivity  E as  e function  of  pL  for  C02 

at  600°K,  based  on  the  empirical  correlations  of  H.  C.  Hottel 
end  R,  B.  Egbert.  Also  shown  ere  the  calculated  contributions 
to  the  total  emissivity  made  by  the  - fundamental. 
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